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We present a novel treatment of resonant massive particles appearing as intermediate states in 
high energy collisions. The approach uses effective field theory methods to treat consistently 
the instability of the intermediate resonant state. As a result gauge invariance is respected in 
every step and calculations can in principle be extended to all orders in perturbation theory, 
the only practical limitation in going to higher orders being the standard difficulties related 
to multi-loop integrals. We believe that the longstanding problem related to the treatment of 
instability of particles is now solved. 



1 The problem with unstable particles 

Physical phenomena studied in many ongoing and future cohider experiments involve the pro- 
duction of heavy unstable particles, such as the top quark t, the vector bosons Z^,W^, the 
Higgs H and, eventually, a range of particles appearing in extensions of the Standard Model. 

One powerful way to determine with high accuracy properties of these particles (i.e. the 
mass M and the width F) is to consider their resonant production. However, if treated as 
stable, unstable particles give rise to non-integrable divergences in cross sections due to the 
internal, divergent propagator i/[q'^ — + ie] (with q the four-momentum of the unstable 
particle). A Dyson-summation of the self-energy n(g, M)^ 

i ^ i ^f U{q,M) Y_ i 

{q-^ - M^) + ie {q'^ - M^) + ie ^ W - My q'^ - - U{q, M) ' ^' 

ensures that the position of the divergence is away from the real axis, since, due to unitarity, 
Im{n(g, M)} 7^ 0. In the neighborhood of = M^, where \Il{q, M)/{q'^ — M^)| > 1, the series 

"Talk given at XXXIX Rencontres de Moriond on Electroweak and Unified Interactions. 



is strictly not convergent. Eq. assumes the vahdity of Perturbation Theory (PT) in the 
whole region and uses analytical continuation in the problematic domain. 

However, eq. resums only a specific class of higher order corrections, namely self-energy 
diagrams, while other loop corrections, such as vertex corrections and box diagrams are neglected 
completely. As a consequence, the resummation procedure jeopardizes all properties valid order 
by order in PT, in particular gauge invariance. Many calculational schemes have been proposed 
in the last years to overcome this longstanding problem (e. g. the fixed width scheme, the 
running width sch eme, the complex mass scheme^, the fermion loop scheme 1^ and the pole 
approximation EEl) . The fundamental problem of most of these approaches is that in order to 
restore gauge invariance physical quantities are modified ad-hoc, introducing a certain degree 
of arbitrariness. Furthermore it is not understood how to improve the accuracy of predictions 
within these frameworks. It turns out that in many cases the size of gauge violating terms is 
numerically small and that various schemes do give quite similar numerical results (though some 
approaches are known to give unreliable predictions at higher center-of-mass energieJ^ . However 
it is important to bear in mind that gauge dependent theoretical predictions have ambiguities 
which can be in principle arbitrarily large if enhanced by suitable kinematical factors or if one 
chooses some pathological gauge. It is therefore mandatory for any precision measurement 
involving resonant unstable particles to develop a framework, relying on first principles only, 
which allows calculations to be performed to the required accuracy while preserving gauge 
invariance. This was the aim of the work presented here. 

2 The model 

To study the conceptual problem of instability it is convenient to consider a simple model which 
allows one to avoid all unnecessary technical difficulties, such as a large number of diagrams, 
integrals containing multiple scales, or the presence of different interactions whose coupling 
constants are numerically different, this would just involve a subtler power counting. 

The toy model that we consider involves therefore only an unstable massive charged scalar 
field, (/>, two massless fermion fields, a charged "electron" ip and a neutral "neutrino" x a-nd 
finally an abelian gauge field, the photon A'^. The scalar decays via Yukawa interaction into an 
electron-neutrino pair. The Lagrangian describing this model is 

+ #^x + y>^xV'-^('/'V)' + ^ct, (2) 

where M and Cct denote the renormalized mass and the counterterm Lagrangian respectively 
and = — igA^ is the covariant derivative. We define Og = /{At[), ay = (yy*)/(47r) (at 
the renormalization scale fi) and assume ~ ~ a, and ax = A/(47r) ~ a^/(47r). 

In the following we will discuss the elements needed to compute the inclusive cross section 

for 

HPi)+e-{p2)^X, (3) 

as a function of = {pi +^2)^ in the resonant region, where S = (g^ — M'^)/q^ ~ a ~ T/M. 
To illustrate the numerics we will use for the MS-mass M = lOOGeV and for the couplings 
ay = ag = 0.1 and a\ = (0.1)^/(47r), evaluated at the renormalization scale /i = M. 

3 The effective theory 

It is well known that in scattering processes involving an unstable particle, radiative corrections 
can be split into factorizable and non-factorizable corrections.!^ The latter are quantum fluctu- 



ations which have a propagation time comparable to the one of the unstable particle r ~ l/F, 
while the propagation time t ~ 1/M ^ r of the factorizable corrections does not allow any quan- 
tum interference with the unstable field. Factorizable corrections are therefore hard loop effects, 
while non-factorizable corrections are described by fluctuations of soft (ps ~ 5M) and collinear 
{pc ~ M,p'^ ~ M5) massless fields or resonant unstable massive particles (g^ — ~ SM). The 
separation between hard and soft /collinear is manifestly gauge invariant and can be performed 
to all orders in PT using an effective field theory approach,'^ rather than a diagrammatic one. 
This classification of radiative corrections can thus be used to construct an effective theory by 
integrating out the hard degrees of freedom which end up in the Wilson coefficients of the ef- 
fective operators, while soft, resonant or collinear fiuctuations remain dynamical modes of the 
effective theory. 

Since all fields have a well-defined scaling in 6, the contributions from different regions can 
be selected with the strategy of regions,!^ which gives an unique prescription of how to expand 
any given loop integral in powers of d. Dimensional regularization ensures than that regions 
with no singularities do not contribute, avoiding the problem of double-counting. Of course, by 
splitting an integral into different regions one might generate spurious singularities. However, 
the infrared singularities of the hard fiuctuations cancels against the ultraviolet singularities of 
the dynamical modes. We choose to subtract them minimally hereby defining a renormalization 
scheme in the effective theory. 

The construction of the effective theory proceeds along the st anda rd way. Let us therefore 
only review the main steps . A ll technical details are given elsewhere.'^ Similarly to Heavy Quark 
Effective Theory (HQET)^2we first redefine the massive field, we extract the big momentum 
component, so that derivatives acting on the new field (p^ produce only soft fluctuations 

= e'^'^^V+cP , d<P, = 0{6) (P.. (4) 

Here V+ projects out the positive frequency states so that is a pure annihilation operator. 
The actual matching procedure is based on requiring that on-shell Greens functions in full theory 
and in the effective theory coincide up to the required accuracy. 

In the resonant region, the effective Lagrangian resulting from eq. ((2)) after integrating out 
the hard modes, is given by a sum of three physically distinct contributions. The propagation 
of the unstable particle gives rise to terms in the effective Lagrangian which closely correspond 
to the well-known HQET Lagrangian 

. / A(^)\ ^ . /(iD [A{i)l2 a(2)\ 
CusET = 2M4 iiv ■ Ds - -^j cP. + 2M4 ( + W -j 

- ^Fsf^uF^ + '^siPs'^s + Xsi^Xs, (5) 

where the subscript s indicates that only the contribution due to soft fields has been included, 
the subscript T indicates the transverse direction with respect to the velocity vector of the 
unstable field. A^*) = O ( a*) are matching coefficients which are defined by 

V A« ^ , (6) 

with s the position of the pole of the propagator, which can be computed from the analytic hard 
part of the self-energy. jChset describes hence the interaction of the unstable particle with soft 
gauge fields and soft fermions. 

The terms in eq. © describing the interaction of energetic fermions in th e in itial and final 
state and gauge fields give rise to the well-known soft-collinear effective theor}^^] 

^SCET = i^c (in-D + il^^,- — i^^^ ^ iP, . (7) 



Here the subscript c indicates contributions from collinear fields along the n_ direction (n^ = 
n\ = 0, n_n+ = 2). Obviously, the SCET Lagrangian for neutrinos can be obtained from eq. 
by replacing covariant derivatives with partial ones. 

Finally, the Lagrangian contains Yukawa vertexes which allow the production and decay of 
the unstable field 

Ant = C y cl)y1pn_Xn+ + C y* (plxn+i^n. + F yy* (V'n_Xn+)(Xn+'0n_) , (8) 

where C = 1 + 0{a) and F are matching coefficients. jCjnt involves new external-collinear fields 
^n-/n+ ^iid Xn_/n+) whosc momenta are given by M /2n_/^ + k, with k = O (6). For instance 
V'n_ is defined as 

The introduction of these fields is useful to distinguish situations were two generic collinear 
fields produce to a state with invariant mass of order M, from situations involving two external 
collinear fields which give rise to states with invariant mass equal to M up to 5M corrections. 

Once the effective Lagrangian has been established, one can compute physical quantities 
using the gauge invariant formalism following from it. In particular we will present here results 
for the inclusive line-shape in our toy model. Since all fields in the effective theory have well- 
defined scaling properties, a power-counting allows to identify the operators that are needed to 
a given order and at which order in the coupling they need to be matched. 

The effective approach is designed to describe the problem in the resonant region, where 
(5 ~ a, therefore at leading order (LO) we resum all terms 0{{a/6)"'). One order beyond, at 
NLO, we include all terms 0{{a/6)^5). NLO terms are expected to correct the LO result by 
~ 5 ~ 10%. In a similar way NNLO corrections resum O 5^) and are expected to be of 

the order ~ 5^ l%. 

At LO the matching amounts simply to resum in the propagator on-shell self-energies. At 
NLO, one needs to match vertex corrections at subleading order and to include sub-leading 
corrections to the propagator. Furthermore dynamical corrections must be also included. 

One of the main results of this work is the possibility to perform a NNLO calculation of 
line-shapes. Even if it is not sensible to struggle with a NNLO calculation in a toy model, it 
is interesting to see what elements are needed to achieve NNLO accuracy. We use the intuitive 
notation where ah denotes an O (a) correction due to a hard loop (and similarly for soft/collinear 
loops). We can thus classify NNLO corrections according to 

• tree amplitude in the effective theory with insertions of LO operators in 5, matched to 
NNLO (al); 

• tree amplitude in the effective theory with insertions of NLO operators in 5, matched to 
NLO 

• tree amplitude in the effective theory with insertions of NNLO operators in 5, matched to 

LO (^2); 

• one-loop amplitude in the effective theory with insertions of LO operators in 5, matched 
to NLO (asO/i, actth); 

• one-loop amplitude in the effective theory with insertions of NLO operators in 5, matched 
to LO {as6, Qc5); 

• two-loop amplitude in the effective theory with insertions of LO operators in 6, matched 
to LO {al, agOc, a^)- 




Figure 1: Left: line shape (in GeV~^) in the MS scheme (solid) and pole scheme (dashed) at LO (upper 
magenta/light grey curves) and NLO (lower blue/dark grey curves) as a function of the center-of-mass energy in 
GeV. Right: The ratio of the NLO to the LO line shape in the MS scheme (solid blue/dark grey curve) and the 
ratio of the NLO MS line shape to a Breit-Wigner fit of the NLO result (dashed magenta/light grey curve). 

This simple classification shows that NNLO calculations of line-shapes are now technically fea- 
sible and that they can be carried out in a systematical and transparent way. We note also that 
due to different kinematics the various contributions are separately gauge invariant. 

We computed the total cross section from the forward scattering amplitude. Fig. ^ (left 
side) shows the LO and NLO line-shape in the MS scheme and pole scheme. The NLO correction 
amounts to 10% at the peak and to up to 30% in the resonant region. The right panel shows the 
ratio of the LO over the NLO line-shape. In the absence of data we performed a fit of the NLO 
line-shape to a Breit-Wigner and we plot the ratio of the NLO result to the Breit-Wigner fit. 
As can be nicely seen from Fig. ^ (right side) the deviation between NLO and the Breit-Wigner 
amounts to up to 15%. Furthermore, the value of the mass fitted from the Breit-Wigner differs 
from the input mass (the true value) by 160 MeV. This shows explicitly that for precision studies 
a proper theoretical prediction should be used, rather than a Breit-Wigner fit, which produce 
sizable shifts in the mass prediction. 



4 Concluding remarks 

The perturbative treatment of unstable particles requires a partial summation of the perturbative 
series, however we believe that the guiding principle was not understood. The breakdown of FT 
is related to the appearance of a second small parameter 6, besides the coupling a. We take the 
attitude that this is the characteristic feature of the problem, so that in a theory that formulates 
this double-expansion correctly, other issues like resummation, gauge invariance and unitarity 
should follow automatically. Since it is a two-scale problem (F <C M) an effective theory is the 
natural framework to formulate such an expansion. The advantages of using effective theory 
methods are 

• calculations are split into well-defined pieces (matching, matrix elements, loop calculations 
in the effective theory . . . ), so that the calculation is efficient and transparent; 

• a power counting scheme in the small parameters {a, 5) allows one to identify terms 
required to achieve a certain accuracy prior to performing the actual calculation; 

• the effective theory provides a set of simpler Feynman rules which allows one to compute 
the minimal set of terms required at the given accuracy. Since one never computes "too 
much", calculations are as simple as possible; 



• calculations can be extended in principle to any order in a, 5, at the price of performing 
complicated, but standard loop integrals; 

• since the expansion has been organized in such a way so as to account for kinematical 
enhancements the PT series in the effective theory converges rapidly; 

• gauge invariance is automatic. 

Despite the simplicity of the model considered (abelian theory, scalar particle), all necessary 
ingredients are provided for the formalism to be applied to any general case. Natural extensions 
concern non-inclusive kinematics, which requires a formalism to expand the real phase-space 
and generally implies that more collinear directions are relevant, and to pair-production of 
unstable particles, in which case the effective Lagrangian will contain two terms similar to 
eq. © describing the propagation of the two particles with different velocity vectors. 
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